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A major  computational  problem  that  arises  in  the  attempt  to  solve 
generalized  network  and  network- related  problems  is  degeneracy.  In 
fact,  using  primal  extreme  point  solution  techniques,  the  number  of 
degenerate  pivots  performed  frequently  ranges  as  high  as  90%  in  large- 
scale  applications.  The  purpose  of  this  paper  is  to  develop  a special  set 
of  structural  and  logical  relationships  for  generalized  network  problems 
that  yields  a new  primal  extreme  point  algorithm  which  avoids  and/or 
exploits  degeneracy.  The  major  mathematical  differences  between  this 
new  algorithm  and  the  simplex  method  are  (1)  each  basis  examined  is 
restricted  to  have  a special  topology,  (2)  the  algorithm  is  finitely 
conve.*'gent  without  reliance  upon  "external”  techniques  such  as  lexicography 
or  perturbation,  and  (3)  a special  screening  criterion  for  nondegenerate 
basis  exchanges  is  available  that  allows  some  of  these  exchanges  to  be 
recognized  prior  to  finding  the  representation  of  an  incoming  arc.  For 
these  reasons,  this  algorithm  has  several  computational  advantages  over 
the  simplex-based  codes  recently  developed  for  solving  generalized  network 
problems. 
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1.  INTRODUCTION 


A generalized  network  (GN)  is  a capacitated  or  uncapacitated  linear  pro- 
gramming (LP)  problem  in  which  the  coefficient  matrix  contains  at  most  two  non- 
zero entries  in  each  column.  Generalized  networks  encompass  a broad  range  of 
important  practical  problems.  Depending  upon  the  non-zero  components  of  the 
coefficient  matrix,  a GN  problem  can  be  a shortest  path  problem,  an  assignment 
problem,  a transportation  problem,  a transshipment  problem,  a generalized  trans- 
portation problem,  or  even  an  optimization  problem  involving  generalized  upper 
bound  (GUB)  conditions  coupled  with  a set  of  disjoint  linear  constraints.  The 
procedures  presented  subsequently  are  general  enough  to  handle  any  of  these 
problems.  To  this  extent,  a higher  degree  of  efficiency  may  be  achieved  for  the 
simpler  problems  by  further  exploiting  their  special  structural  characteristics. 
The  more  complex  variants  of  the  GN  class  have  received  the  attention  of  several 
studies  [2,  3,  4,  11].  The  prevalence  of  these  problems  in  practical  settings 
[6,  7,  8,  9,  10,  12,  17,  20,  21]  has  created  a need  for  new  theoretical  results 
that  can  provide  efficient  solution  methods. 

A major  computational  problem  that  arises  in  the  attempt  to  solve  GN  and 
GN-related  problems  is  degeneracy.  In  fact,  using  primal  extreme  point  solution 
techniques,  which  are  currently  the  most  efficient  known  for  these  problems,  the 
number  of  degenerate  pivot  steps  is  frequently  as  high  as  90%  in  large-scale 
applications.  The  purpose  of  this  paper  is  to  develop  a special  set  of  structural 
and  logical  relationships  for  GN  problems  that  yields  a new  primal  extreme  point 
algorithm  which  avoids  and/or  exploits  degeneracy.  The  conceptual  foundations 
underlying  this  algorithm  may  be  viewed  as  an  extension  of  those  underlying  the 
recently  developed  algorithms  of  [2,  3,  4,  11]  for  solving  assignment,  transporta- 
tion, and  transshipment  problems.  However,  the  more  complex  structure  of  gen- 
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eralized  networks  requires  several  theoretical  departures  and  advances.  One  of  the 
principal  features  of  the  new  algorithm,  shared  In  common  with  Its  earlier  cousins 
for  less  general  problems.  Is  a strong  form  of  convergence  that  limits  the  number 
of  degenerate  steps  In  a far  more  powerful  way  than  achieved  by  "lexicographic 
Improvement"  (as  used,  for  example.  In  customary  LP  perturbation  schemes). 

Each  basis  examined  by  this  algorithm  Is  restricted  to  have  a special 
topology.  We  show  that  If  a GN  problem  has  an  optimal  solution,  then  an  optimal 
solution  can  be  found  by  considering  only  bases  of  this  type.  The  major  mathe- 
matical differences  between  this  new  algorithm  and  the  simplex  method  are  (1)  the 
rules  of  the  algorithm  automatically  (without  search)  ensure  that  all  bases  have 
the  special  topological  structure,  and  bypass  all  other  bases  normally  given  con- 
sideration by  the  simplex  method;  (2)  the  algorithm  Is  finitely  convergent  without 
reliance  upon  "external"  techniques  (such  as  lexicography  or  perturbation);  and 
(3)  a special  screening  criterion  for  nondegenerate  basis  exchanges  Is  available 
that  allows  some  of  these  exchanges  to  be  recognized  prior  to  finding  the  repre- 
sentation of  an  Incoming  arc.  For  these  reasons,  this  algorithm  has  several 
computational  advantages  over  the  simplex-based  codes  recently  developed  for 
solving  GN  problems. 


— 1 
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2.  PROBLEM  STATEMENT  ! 

I 

1 

The  generalized  network  problem  may  be  defined  as  follows:  j 

Primal  i 


Minimize 
Subject  to 

T 

C X 

(1) 

AX  = b 

(2) 

0<X<U 

(3) 

Maximize 

T,^  T 

TT  b - Y U 

(M 

Subject  to 

T T 

TT  A - 

(5) 

Ti  - unrestricted 

(6) 

o 

A 1 

(7) 

where  U > 0 and  A is  an  m x n matrix  in  which  each  column  contains  at  most  two 
j nonzero  coefficients  of  opposite  sign.  We  will  assume  that  one  coefficient  is  -1 

and  the  other  any  positive  real  number.  When  a column  has  only  one  nonzero  coef- 
ficient, we  will  assume  that  this  coefficient  is  either  1 or  -1.  A variable 
associated  with  such  a column  is  called  a slack  variable  if  its  nonzero  coefficient 
is  1,  and  is  called  a surplus  variable  if  its  nonzero  coefficient  is  -1. 

The  most  efficient  procedures  for  solving  GN  problems  [13,  14,  19,  22,  24] 
are  based  on  viewing  the  problem  in  a graphical  context.  These  procedures  are 
adaptations  of  the  simplex  algorithm  [1,  7,  12]  in  which  the  A matrix  and  the 
basis  matrix  are  stored  as  graphs  using  computer  list  structures.  The  use  of  such 
structures  reduces  both  the  amount  of  work  needed  to  perform  basic  simplex  opera- 
tions and  the  amount  of  computer  memory  required  to  store  essential  problem  data. 
Computation  and  memory  are  further  reduced  if  A does  not  have  full  row  rank,  since 
such  GN  problems  are  equivalent  to  a disjoint  set  of  transshipment  problems  [15]. 
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In  addition,  the  graphs  contain  only  the  nonzero  matrix  components  which  allow 
special  graph  labeling  procedures  [18,  22,  24]  to  eliminate  checking  and  un- 
necessary arithmetic  operations  on  zero  elements. 

The  development  of  the  algorithm  in  this  paper  demonstrates  an  additional 
advantage  of  representing  and  solving  GN  problems  graphically.  We  will  show 
that  the  graphical  representation  allows  one  to  fully  characterize  the  nonzero 
elements  of  the  representation  of  any  nonbaslc  vector  and  the  signs  of  these 
elements.  We  will  also  show  that  the  representation  allows  one  to  explicitly 
characterize  which  dual  values  change  and  how  these  values  are  altered  after 
performing  a basis  exchange  step.  These  mathematical  characterizations  and 
visualizations  provide  the  cornerstones  of  the  development  and  proof  of  our 
algorithm. 

The  essential  concepts  and  definitions  of  elementary  graph  theory  that  will 
be  used  in  our  development  are  presented  and  related  to  the  GN  problem  in  the 
next  section. 

3.  BACKGROUND 

A directed  graph  G(N,  E)  is  a finite  set  of  nodes  N and  arcs  E connecting 
the  nodes.  Each  element  of  E is  an  ordered  pair  (i,  j)  of  distinct  elements  of 
N which  represents  an  arc  from  the  node  i to  node  j.  If  the  arc  set  E is  ex- 
panded to  contain  arcs  which  have  both  endpoints  incident  on  the  same  node,  or 
multiple  arcs  connecting  the  same  two  nodes,  then  G(N,  E)  is  called  a general 
graph. 

The  GN  problem  defines  a general  graph  as  follows.  Each  row  of  A corres- 
ponds to  a node  and  each  column  corresponds  to  an  arc  of  the  graph.  The  nonzero 
entries  in  a column  will  be  called  the  multipliers  of  the  corresponding  arc.  Asso- 
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elated  with  each  arc  is  a variable,  an  upper  bound,  an  objective  function  co- 
efficient, and  the  nonzero  multipliers.  The  variable  (i.e.,  the  component  of  X) 
associated  with  the  arc  is  called  the  flow  on  the  arc.  (Or,  frequently,  the  value 
of  this  variable  is  called  the  flow.)  The  row  positions  of  the  multipliers  in 
the  column  that  corresponds  to  the  arc  identify  the  nodes  on  which  the  arc  is 
incident.  The  arc  is  directed  from  the  node  associated  with  the  -1  multiplier  to 
the  node  associated  with  the  positive  multiplier. 

An  arc  directed  from  node  i to  node  j will  be  denoted  (i,  j).  The  positive 
multiplier  associated  with  arc  (i,  j)  will  be  denoted  (Since  there  may  be 

more  than  one  arc  (i,  j)  from  i to  j , technically  a third  index  should  be  used; 
however,  for  notatlonal  convenience,  the  third  subscript  will  be  omitted.  The 
method  as  subsequently  described  provides  an  organization  by  which  multiple  arcs 
with  unique  multipliers,  costs  and  capacities,  are  readily  accommodated.)  If  a 
column  has  only  one  multiplier  (nonzero  entry)  both  endpoints  of  the  arc  are 
incident  on  the  same  node  and  the  arc  is  said  to  form  a loop.  Such  an  arc,  with 
endpoints  incident  on  the  same  node  i,  is  customarily  denoted  (i,  i)  and  its 

multiplier  is  denoted  a... 

11 

The  right-hand  side  vector  b for  the  GN  problem  associates  a node  requirement 
bj^  (the  component  of  b)  with  node  k of  the  network.  (Each  unit  of  flow  on  an 
arc  (i,  j)  therefore  "contributes"  -1  and  , respectively,  to  the  node  require- 
ments b.  and  b..)  In  addition,  each  node  k has  an  associated  dual  variable  tt  , 

1 J k 

called  its  node  potential. 


4.  BASIS 


4.1  Basis  Structure 


Using  the  standard  bounded  variable  simplex  algorithm,  a basis  B for  the  GN 
problem  is  a full  row  rank  matrix  composed  of  a linearly  independent  set  of 


column  vectors  of  A (or  possibly  of  A augmented  by  unit  vectors  corresponding  to 
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artificial  variables).  The  variables  associated  with  the  columns  of  B are  considered 
to  be  basic  variables  and  all  others  are  nonbasic  variables.  A basic  solution 
consists  of  assigning  each  nonbasic  variable  a value  equal  to  its  lower  or  upper 
bound.  Thereupon,  a unique  value  results  for  each  basic  variable  to  satisfy 
equation  (2). 

A basis  for  a GN  problem  may  be  viewed  and  stored  as  a graph  which  contains 
only  the  nonzero  components  of  the  matrix  B.  The  characteristics  of  the  graph  cor- 
responding to  B have  been  fully  outlined  in  the  literature  [1,  12,  23]  when  the 
graph  of  A is  bipartite.  Several  authors  [18,  22,  24]  have  further  noted  that  the 
structure  of  the  graph  corresponding  to  B is  essentially  the  same  (except  for  arc 
orientations)  when  the  graph  of  A is  not  bipartite.  For  completeness  we  will 
rigorously  characterize  the  structure  for  the  non-bipartite  case. 

Clearly,  by  the  correspondences  already  indicated,  any  basis  for  the  GN 
problem  may  be  viewed  as  a graph.  The  basis  graph  contains  all  of  the  nodes  but 
only  a subset  of  the  arcs  of  the  graph  of  A,  and,  hence  is  called  a spanning 
graph.  A spanning  graph  consists  of  one  or  more  connected  subgraphs,  where  a 
connected  graph  (hence  subgraph)  is  a graph  in  which  a path  exists  between  any 
two  distinct  nodes. 

th 

Let  B^  denote  the  i connected  subgraph  of  B.  Then  B can  be  represented  as 
a block  diagonal  matrix  as  follows: 


B 

n 


nilr’i  -I  • nl  v 
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By  the  non-singularity  of  B,  each  is  a square  n^  x n^  non-singular  matrix  and 
consequently  is  a connected  graph  consisting  of  n.  nodes  and  n,  arcs.  Therefore, 
it  follows  that  the  subgraph  corresponding  to  each  B^  is  a quasi-tree  [18],  i.e., 
a simple  tree  to  which  a single  arc  has  been  added.  The  additional  arc  creates 
exactly  one  loop  (a  circular  path)  in  the  quasi-tree.  As  noted  earlier,  an  arc 

/ 

having  both  endpoints  incident  on  the  same  node  (i.e.,  a slack  or  surplus  variable) 
forms  a loop  henceforth  referred  to  as  a self-loop . The  basis  structure  may  thus 
be  simply  characterized  as  follows. 

/ 

Remark  1 : By  arranging  rows  and  columns,  any  basis  B for  a GN  problem  is  a 

block  diagonal  matrix  whose  blocks  are  triangular  except  for  rows  and  columns 
corresponding  to  loop  arcs  (i.e.,  arcs  belonging  to  loops  asso^ated  with  the 
blocks).  A block  is  completely  triangular  if  its  loop  is  a/self-loop. 

Proof : The  proof  follows  immediately  from  the  fact  tha^r  a basis  B for  a GN 

problem  consists  of  one  or  more  disjoint  quasi-trees. 

An  efficient  method  called  the  extended  a^mented  predecessor  index  (EAPI) 
method  [18]  has  been  developed  for  storyjg  and  updating  the  basis  graph  of  the  GN 
problem.  One  of  the  advantages  of Jfne  EAPI  method  is  its  ability  to  find  the 
representation  of  an  entering/Wctor  and  to  determine  updated  values  for  dual 
variables  by  tracing  app^priate  segments  of  the  quasi-trees  of  the  basis  graph. 
Therefore,  no  expT^lt  basis  inverse  is  needed,  and  tiie  usual  arithmetic  operations 
required  to  update  the  inverse  are  eliminated.  (The  accumulition  of  roundoff  errors 
is  likewise  substantially  reduced  or  eliminated). 

The  EAPI  method  also  plays  an  important  role  in  the  statement  and  develop- 
ment of  the  new  extreme  point  algorithm  presented  in  this  paper.  Its  special  up- 
dating rules  can  be  conveniently  expressed  in  terms  of  a set  of  operations  applied 
to  Ellis  Johnson's  triple  label  representation  for  recording  and  tracing  a tree, 
adapted  to  accommodate  structural  peculiarities  of  a quasi-tree.  The  approach 
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assigns  three  labels  to  each  node  of  a quasi-tree:  a predecessor,  a successor,  and 

a brother.  This  triple-label  scheme  is  first  applied  to  the  portion  of  each  quasi- 
tree that  results  by  suppressing  all  loop  arcs  of  the  quasi-tree  so  that  the  labeled 
portion  is  a disjoint  set  of  rooted  trees.  The  root  of  each  tree  consists  of  the 
node  that  lies  on  the  loop;  that  is,  the  arcs  are  oriented  by  the  predecessor  labels 
so  that  the  unique  path  from  any  node  to  the  root  node  of  the  tree  is  a directed 
path.  (In  this  manner  the  predecessor  labels  uniformly  orient  the  tree  to  create 
an  arborescence . ) Notationally , we  let  p(i)  denote  the  predecessor  of  node  1. 

The  triple-label  representation  may  be  viewed  as  inducing  an  "ancestry  rela- 
tionship" on  a tree,  each  node  carrying  three  node  indexes,  which  name  the  father 
(predecessor),  the  eldest  son  ("first"  successor),  and  next  younger  brother  of  the 
given  node.  A node  is  taken  to  be  the  father  of  all  its  immediate  successors,  these 
latter  constituting  a set  of  brothers,  arbitrarily  sequenced  from  eldest  to  youngest. 
Thus,  the  root  node  is  the  common  ancestor  of  all  nodes  in  the  tree.  Nodes  at  the 
extremities  of  the  tree  have  no  immediate  successors  and  the  "last"  of  a set  of 
successors  of  a given  node  has  no  younger  brother.  In  these  extreme  cases,  a 
"dummy"  name  which  communicates  their  nonexistence  is  used. 

The  EAPI  method  connects  these  trees  to  the  remainder  of  the  quasi-trees  in 
which  they  lie  by  reinstating  the  loop  arcs.  Then  the  root  of  each  tree  is  as- 
signed a predecessor  which  orients  these  loop  arcs  uniformly  clockwise  or  counter- 
clockwise. Each  node  on  the  loop  thus  becomes  its  own  ancestor.  (A  self-loop  may 
simply  be  assigned  the  orienration  it  receives  in  the  graph  of  A.)  Hence,  each 
loop  node  has  an  equal  status  as  an  ancestor  of  all  nodes  ir  the  quasi-tree.  This 
is  called  the  rooted  loop  orientation  of  a quasi-tree  [18]. 

The  basis  graph  illustrated  in  Figure  1 consists  of  two  quasi-trees.  The  arcs 
are  drawn  according  to  their  actual  direction  in  the  graph  of  A rather  than  to  the 
orientation  imparted  by  the  predecessor  indexing.  The  basis  matrix  B associated 


Figure  1 - Basis  Graph  of  A 


Figure  2 - EAPI  Basis  Graph 
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Table  I - EAPI  Representation  of  a Basis 
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with  this  graph  consists  of  a linearly  independent  set  of  column  vectors  from  A 
(possibly  augmented  by  "artificial"  unit  vectors).  The  same  basis  graph  repre- 
sented in  EAPl  form  where  each  quasi-tree  possesses  a rooted  loop  orient.ation  is 
illustrated  in  Figure  2. 

To  formulate  and  derive  the  results  of  this  paper  in  the  simplest  fashion,  it 

is  convenient  to  Introduce  the  notion  of  a mirror  arc.  (This  notion,  in  less 

general  form,  is  Introduced  in  the  pure  network  setting  in  [5]).  The  mirror  of  an 

arc  (i,  j),  with  cost  c..,  multiplier  a , and  flow  x satisfying  u..  > x >0, 

iJ  ij  ij  ij  - iJ  - 

is  an  oppositely  directed  arc  (i,  i)  with  cost  -c  /a  , multiplier  1/a  , and 

ij  1.1  ij 

flow  X..  satisfying  -u . . a..  < x..  <0.  The  mirror  arc  is  clearly  the  result  of  a 

simple  linear  transformation  by  which  x . = -x  ’a  . Viewed  alternatively,  the 

Ji  iJ  iJ 

mirror  arc  arises  by  a scaling  operation  that  divides  the  column  of  A associated 

with  the  variable  x..  by  the  quantity  -a.  . Thus  the  multiplier  a.,  that  lies  in 
iJ  ij  ij 

row  j is  changed  to  -1,  indicating  that  the  mirror  arc  will  be  directed  out  of, 
rather  than  into,  node  ].  Similarly,  the  -1  in  row  i becomes  1/a^^,  identifying 
this  to  be  the  multiplier  of  the  new  arc.  The  scaling  is  applied  correspondingly 
to  the  appropriate  components  of  the  vectors  c and  u.  Note  that  there  is  a con- 
venient reciprocal  relationship  between  an  arc  and  its  mirror.  That  is,  if  the 
same  type  of  scaling  process  is  applied  a second  time,  we  discover  that  the  mirror 
of  the  mirror  is  the  original  arc,  or  in  other  words,  each  arc  is  the  mirror  of  the 

other.  (By  further  employing  the  linear  translation  x.  = -x.,/a..  + u , both 

1,1  U iJ  11 

arcs  would  be  constrained  to  nonnegative  flows.  We  do  not  bother  with  the  trans- 
lation since  it  requires  an  extra  step  and  is  artificial  for  arcs  with  infinite 
capacit ies . ) 

As  a consequence  of  these  remarks,  an  arc  may  be  replaced  (conceptually  or 
literally)  at  any  time  by  its  mirror,  provided  the  appropriate  relationships  are 


'•»  i • - I ^ 
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observed.  Tluis , for  example,  a flow  increase  of  1 unit  on  arc  (i,  j)  corresponds 
to  a flow  decrease  of  a^.  units  on  its  mirror  arc  ( j , i),  and  if  the  flow  on  arc 

(i,  j)  is  y units  below  its  upper  bound,  then  the  flow  on  the  mirror  arc  ( j , i)  is 

ya^^  units  above  its  lower  bound.  (As  noted  earlier,  we  allow  the  existence  of 
multiple  arcs  between  the  same  two  endpoints,  but  continue  to  refer  to  a given 
arc  from  i to  j as  simply  (i,  j),  suppressing  a third  subscript  for  notational 
convenience.  [The  word  "given"  implicitly  provides  the  third  subscript.]) 

The  notion  of  a mirror  arc  can  be  applied  to  slack  and  surplus  arcs  by  means 
of  a refinement  of  the  self-loop  concept.  A self-loop  at  a node  i is  conceived 
as  joining  node  i to  an  implicit  copy  of  i,  denoted  i#.  A surplus  arc  at  node  i 
(the  self-loop  that  has  a -1  in  the  i^^  position  of  its  column)  will  be  designated 

as  (i,  i#) , and  a slack  arc  at  node  i (the  self-loop  that  has  a 1 in  the  i*"'^  posi- 

tion of  its  column)  will  be  designated  as  (i#,  i).  The  multiplier  for  both  a 
slack  arc  and  surplus  arc  is  taken  to  be  1.  Consequently,  as  a result  of  this 
refinement  and  the  previous  definitions,  the  mirror  of  an  arc  (i,  i#) , with  cost 

c.,„,  and  flow  x..,,  satisfying  u . . > x..,,  > 0,  is  the  arc  (i#,  i)  with  cost 
11#  11#  11#  - 11#  - 

c.j, . = -c  . . and  flow  x . . satisfying  -u . , , < x . . < 0.  (The  mirror  of  an  arc 
(i#,  i)  is  an  arc  (i,  i#)  characterized  in  a like  fashion.) 

The  advantage  of  introducing  the  notion  of  a mirror  arc  is  twofold.  The  mirror 
arc  concept  makes  it  possible  to  ensure:  (1)  all  nonbasic  arcs  have  flows  equal 

to  their  lower  bounds,  and  (2)  all  basic  arcs  have  the  same  direction  as  imparted 
by  the  orientation  of  the  predecessor  indexing.  A nonbasic  arc  whose  flow  equals 
its  upper  bound  or  a basic  arc  whose  direction  is  contrary  to  the  predecessor 
orientation  may  simply  be  replaced  by  its  mirror  arc,  thereby  enabling  both  (1) 
and  (2)  to  hold.  It  may  be  noted  that  a slack  arc  (i#,  i)  in  the  basis  automati- 
cally is  directed  in  accordance  with  the  orientation  of  the  predecessor  indexing, 


while  a surplus  arc  (i,  i#)  is  directed  contrary  to  this  orientation,  requiring  the 
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latter  to  be  replaced  by  its  mirror.  (This  follows  from  the  natural  convention 
that  the  condition  p(i)  = i for  a self-loop  implicitly  represents  the  condition 
p(i)  = i//.)  This  replacement  of  arcs  by  their  mirrors  to  give  them  a desired 
orientation  or  a flow  equal  to  a lower  bound,  has  important  conceptual  advantages, 
as  will  become  clear  in  the  subsequent  development. 

Hereafter,  we  will  continue  to  use  the  term  mirror  of  an  arc  to  refer  to  an 
arc  obtained  from  any  given  arc  (original  or  otherwise)  by  the  mirror  operation, 
but  will  restrict  the  term  mirror  arc  to  refer  to  the  mirror  of  an  original  arc. 

From  the  foregoing,  the  basis  arc  (p(i),  i)  , directed  to  node  i from  its 
predecessor  node  p(i),  is  an  original  arc  if  the  associated  basic  column  of  A has 
a positive  coefficient  in  position  i and  is  a mirror  arc,  otherwise.  The  scaled 
matrix  that  results  from  a basis  B in  this  fashion  will  be  denoted  B'.  (That  is, 
each  column  of  B'  gives  the  appropriate  representation  of  the  corresponding  basis 
arc . ) 

Some  notational  simplification  is  now  possible.  Since  each  node  i in  the 
basis  g'aph  associated  with  B'  has  a unique  predecessor,  the  unique  multiplier  for 
arc  (p(i),  i)  (which  may  be  either  an  original  arc  or  a mirror  arc)  can  be  asso- 
ciated with  node  1 and  will,  therefore,  be  den)ted  by  a^.  In  addition,  the  unique 
variable  corresponding  to  each  arc  (p(i),  i)  can  be  associated  with  node  i by  re- 
ordering the  columns  of  B'  so  that  the  i*"*^  column  vector  in  B'  is  associated  with 
arc  (p(i),  i) . Thus  the  basic  variable  associated  with  arc  (p(i),  i)  in  the 
. basis  graph  will  be  denoted  by  x^.  This  is  illustrated  in  Figure  2. 

Due  to  this  scaling  and  reordering,  the  information  required  to  represent  a 
basis  in  EAPI  form  can  be  organized  in  an  efficient  manner  as  shown  in  Table  1. 
Thus,  in  Table  I,  the  array  NODE  contains  node  name  and  would  not  be  stored. 

Arrays  FRED,  SUC,  and  BRO  contain  the  predecessor,  successor,  and  brother  labels, 
respectively,  associated  with  each  node  name.  A zero  label  is  used  to  indicate  the 


j 


nonexistence  of  a label.  Associated  with  each  node  i in  the  arrays  FLOW  and 


MULT  is  the  flow  and  multiplier,  respectively,  on  the  arc  (p(i),  i)  from  the 
predecessor  of  i to  i.  Also  associated  with  each  node  i in  the  MIRROR  array  is  a 
value  of  1 if  the  arc  (p(i),  i)  is  a mirror  arc  and  0 if  it  is  an  original  arc. 

This  organization  is  in  effect  a compact  storage  scheme  for  B since  B can  be 
created  using  the  arrays  PRED,  MULT,  and  MIRROR. 

The  conventions  described  not  only  facilitate  an  effective  computer  programming 
implementation,  but  greatly  simplify  the  notation  required  to  express  and  prove  a 
number  of  our  main  results.  It  is  important  to  keep  in  mind  that  results  derived 
graphically  (in  terms  of  arcs  of  the  EAPI  basis  graph)  have  a direct  algebraic 
analog  (in  terms  of  the  basis  matrix  B'  associated  with  this  graph).  Likewise, 
results  derived  algebraically  by  reference  to  B'  translate  immediately  to  the  EAPI 
basis  graph.  We  will  find  it  convenient,  at  various  points,  to  alternate  the  type 
of  derivation  employed,  though  we  will  continue  to  rely  on  the  basis  graph  as  the 
chief  vehicle  for  motivating  and  elucidating  our  results. 


4 . 2 Exchange  Steps  and  Basis  Representations 

The  procedures  of  an  adjacent  extreme  point  algorithm  identify  an  incoming 
arc  (p,  q)  and  an  outgoing  arc  (r,  s),  respectively,  from  among  the  nonbasic  and 
basic  arcs.  It  should  be  pointed  out  that  (p,  q)  can  refer  to  a slack  arc  (p  = q#) 
or  a surplus  arc  (q  = p#) , and  (r,  s)  can  refer  to  a slack  arc.  (The  outgoing  arc 
cannot  refer  to  a surplus  arc  while  it  is  in  the  basis,  in  order  to  retain  proper 
orientation,  as  observed  previously.)  The  addition  of  (p,  q)  and  the  removal  of 
(r,  s)  from  the  current  basis  produces  a changed  set  of  quasi-trees,  and  this  opera- 
tion constitutes  a fundamental  step  or  "basis  exchange"  step  of  linear  programming 
methods . 


r 
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Basis  Representat Ion  of  an  Incoming  Arc  j 

In  any  quasl-tree  possessing  a rooted  loop  orientation.  It  Is  clear  that  a 
sequential  trace  of  predecessors  of  a given  node  generates  a predecessor  path  which 

1 

contains  all  arcs  on  the  loop.  The  arcs  of  a predecessor  path  are  encountered  In 
the  sequence  determined  bv  tracing  from  descendant  to  ancestor,  rather  than  from 
ancestor  to  descendant  (though  the  latter  trace  accords  with  the  orientation  of 
the  basis).  For  simplicity,  we  shall  suppose  that  such  a path  Is  simple;  l.e., 
not  traced  beyond  necessity.  Hence,  as  soon  as  a node  on  this  path  Is  encountered 
a second  time,  the  trace  Is  stopped  (which  applies  also  to  the  situation  p(l)  = 1, 

Implicitly  representing  p(l)  = 1//).  The  path  Itself  therefore  duplicates  no  arcs. 
Henceforth,  will  denote  the  set  of  arcs  In  the  predecessor  path  of  node  i. 

We  will  follow  the  convention  that  the  predecessor  path  for  a node  which  Is 

an  Implicit  copy  of  another  Is  empty.  Thus,  In  particular.  If  an  arc  (1,  j ) Is  a 
self-loop  then  = 0 If  (1,  j)  represents  a slack  arc  (1  = j//)  and  P^  = 0 If 

(1,  j)  represents  a surplus  arc  (j  = 1//).  (Note  that  P^  Is  never  empty  for  any 

node  1 that  Is  not  an  Implicit  copy  of  another,  since  the  predecessor  path 

must  always  end  In  a loop  that  contains  at  least  one  arc.) 

An  Important  aspect  of  a basis  exchange  step  Is  to  Identify  the  set  of  arcs 
on  which  flows  will  change,  or  more  precisely,  the  set  of  basic  arcs  which  provide 

a linear  representation  of  the  nonbaslc  Incoming  vector.  The  column  vector  cor- 

responding to  the  Incoming  arc  (p,  q)  when  (p,  q)  does  not  correspond  to  a self- 
loop Is  of  the  form  -E^  + a E*^,  where  E^  denotes  the  1^^  column  of  the  Identity 

pq 

matrix.  If  the  Incoming  arc  (p,  q)  corresponds  to  a self-loop,  the  column  vector 
corresponding  to  this  arc  Is  of  the  form  E*^  (If  p = q#)  or  -E^  (If  q = p#) . It 
Is,  therefore,  convenient  to  first  characterize  the  representation  of  any  unit 
vector  E^  with  respect  to  the  arcs  In  the  basis  graph. 

t 

I 

A 


II  wid  I HI  '\k  1^ 
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Remark  2 : The  basic  arcs  (vectors)  that  have  a nonzero  coefficient  in  the  repre- 

sentation of  a unit  vector  are  contained  in  P . 

1 

Instead  of  proving  this  remark  in  its  present  form,  we  will  state  and  construc- 
tively prove  a more  encompassing  result  that  illustrates  several  useful  principles. 
To  lay  the  foundation  for  this  result,  we  suppose  the  cardinality  of  is  m. 

Then  P.  is  a set  of  m arcs  which  connect  m distinct  nodes.  We  renumber  these 
1 

nodes  so  that  i = 1 and  the  predecessor  of  node  k is  k+1  for  k = 1,  . . . , m-1. 
(This  corresponds  to  reordering  the  rows  of  the  matrix  B'  associated  with  the 
basis  graph.)  For  the  renumbered  nodes  let  node  w (w  > 1)  denote  the  first  loop 
node  encountered  in  a sequential  trace  of  the  predecessors  from  node  i.  (Note  that 
m > w and  the  predecessor  of  m is  the  duplicate  node  encountered  in  the  predecessor 
path  from  node  1 and  thus  is  node  w.)  Figure  3 graphically  illustrates  P^  subject 
to  this  node  numbering.  We  also  suppose  that  the  arcs  in  P^  are  numbered  cor- 
respondingly, so  that  the  k*"^  arc  is  (p(k),  k)  for  k = 1,  . . . , m.  Using  the 
notation  established  in  the  previous  section,  the  multipliers  and  columns  of  B' 

associated  with  these  arcs  are  a, , . . . , a and  B'  , . . . , B'  , respectively, 

1 ml  m 

t h 

where  B',  denotes  the  k column  of  B'. 
k 

Remark  2 asserts  that  there  exists  a solution  to  the  equation 

B'.Y,  + B’-Y^  + . . . . B'  Y = (8) 

11  2 2 ram 

For  our  purposes,  however,  we  will  go  beyond  establishing  the  existence  of  such 
a solution  and  identify  its  form  exactly. 

For  a self-loop,  where  m = w,  equation  (8)  represents  the  following 
triangular  system: 


Where  = 1 (since  the  multiplier  of  a self-loop  is  always  1). 

For  other  than  a self-loop,  equation  (8)  represents  the  following  which  is 

identical  to  the  preceding  in  the  first  w-1  components,  but  which  differs  over 
I 

the  "loop  portion:" 


Let  = the  set  of  the  first  k arcs  encountered  by  a trace  of  the  predecessor 
path  from  node  1 (Sj^  = {(p(i),  1),  i = 1 . . . , k})  and  let  L = the  set  of  arcs  on 
the  loop  of  the  predecessor  path  from  node  1 (L  = {(p(i),  i) , i = w,  . . . , m}). 

In  general,  for  any  set  of  arcs  S,  we  define 

G (S)  = l/(the  product  of  the  multipliers  of  the  arcs  in  S). 
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Thus,  in  particular 


e(Sj^)  = l/a^a^  • • . 


0(L)  = 1/a 


The  quantity  0(L)  is  called  the  loop  factor  of  the  loop  L.  (By  convention,  we 
define  0(0)  = 1,  where  0 is  the  empty  set.) 

Remark  3:  The  solution  to  (9)  is  given  by 
= 0(Sj^)  for  k = 1,  . . . , w 
and  the  solution  to  (10)  is  given  by 

y^^  = k = 1,  . . . , w-1 

yj^  = 0(Sj^)/(l  -0(L))  for  k = w,  . . . , m 

Proof : The  solution  to  the  triangular  system  (9)  is  clearlj 
= l/a^^  = 0(S^) 

^k  = Vl^\  = for  k = 2,  . . . , w 

Similarly,  the  solution  to  the  triangular  portion  of  (10)  is  the  same  as  for  (9) 

for  k = 1,  . . . , w-1,  and  it  is  only  required  to  show  y^^  = 0(Sj^)/(l  -0(L))  for 

k = w,  . . . , m.  The  first  equation  of  the  loop  portion  of  (10)  yields 

y = (y  1 y )/a 

w w-1  m w 

and  subsequent  equations  yield 

^k  = Vl^\  = ^Vl  for  k = w+1,  . . . , m 

Hence  in  particular, 

y = (y  1 y )/a  . . . a = (0(S  ,)  + y )’0(L))  = 

m w-1  m w m w-i  m 

0(S  ,)  0(L)  + y 0(L)  = 0(S  ) + y 0(L) 

w-1  m m m 

Solving  for  y gives 
m 

y = 0(S  )/(l-0(I,)) 

m IT) 


J 
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The  remark  then  follows  from  the  fact  that  0(S,  ) = 0(S,  ,)/a,  and  y, 

k k-1  k k 


= Vi/\ 


In  terms  of  the  EAPI  graph,  the  solution  Identified  in  Remark  3 can  be 
viewed  as  the  determination  of  flows  on  the  arcs  of  to  satisfy  a requirement  of 
1 at  node  i.  The  values  k = 1,  . . . , m,  identify  these  arc  flows.  Further, 
the  way  in  which  the  value  of  y^^  is  obtained  from  the  value  of  y^^  over  the 
triangular  portion  of  the  system  can  be  viewed  as  the  transmission  of  the  require- 
ment at  node  i to  successive  nodes  of  as  a consequence  of  assigning  flow  values 
along  the  path.  Thus,  in  particular,  the  graphical  interpretation  applied  to  the 

triangular  portion  of  the  system  can  be  expressed  in  general  as  follows: 

s t 

The  assignment  of  a flow  value  of  to  the  (k-1)  arc  of  the  path  transmits 

a requirement  equal  to  to  the  k*"^  node  of  the  path.  In  turn,  this  require- 

ment of  y^^  ^ at  the  k*"^  node  of  the  path  determines  the  unique  flow  value  of 

yj^  = ^k-l^^k  ~ ^^^1  ■ ■ ‘ ^k  ~ path.  The 

solution  for  the  loop  portion  of  the  system  (in  the  non-triangular  case)  can  be 

viewed  similarly.  The  assignment  of  flows  to  the  non-loop  portion  of  the  system 

transmits  a requirement  of  y = 6(S  , ) to  node  w (equation  w)  and  this  must  be 

satisfied  by  assigning  flows  to  the  arcs  (values  to  the  variables)  associated  with 

the  loop.  Specifically,  as  shown  in  the  proof  of  Remark  3,  this  yields  a flow  of 

y = 0(S  )/(l  -0(L))  on  the  first  arc  of  the  loop,  whereupon  the  remaining  flows 
w w 

and  transmitted  requirements  are  identified  exactly  as  in  the  triangular  case, 

i-e.,  yj^  = for  k = w+1,  . . . , m. 

The  flows  thus  determined  provide  the  representation  for  the  unit  vector 

relative  to  B'.  These  must  be  re-scaled  to  obtain  the  representation  relative  to 

t h 

the  basis  matrix  B.  Thus  the  k element  of  the  re-scaled  solution  equals 
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if  the  k arc  is  an  original  arc,  and  equals 
. _ , , th 

It  the  k arc  is  a mirror  arc 
for  k = 1,  . . . , m. 

Remark  4:  The  predecessor  paths  P and  P , for  two  distinct  nodes  p and  q. 

- p q 

contain  all  basic  arcs  with  a nonzero  coefficient  in  the  basis  representation  of 
arc  (p,  q),  where  arc  (p,  q)  represents  any  one  of  the  possible  multiple  arcs 
connecting  node  p to  node  q. 

Proof : This  remark  follows  directly  from  Remark  2 since  finding  the  representa- 

tion of  arc  (p,  q)  is  equivalent  to  finding  -E^  + a if  (p,  q)  does  not  corres- 

pq 

pond  to  a self-loop.  Otherwise,  the  representation  is  equivalent  to  finding  E^  if 


p = q#  and  to  finding  -E^  if  q = pit. 

The  foregoing  procedures  make  it  possible  to  determine  the  representation 


of  a nonbasic  arc  (p,  q)  graphically  as  follows;  Find  appropriate  flows  on  the 
arcs  in  Pp  (if  p ^ qit)  which  satisfy  a node  requirement  at  node  p of  -1.  Likewise, 
find  appropriate  flows  on  arcs  in  P^  (if  q ^ pit)  which  satisfy  a node  requirement 
at  node  q of  Assign  zero  flows  to  all  other  basic  arcs.  Add  the  resulting 

sets  of  flows  to  yield  the  coefficients  which  provide  the  representation  of  a non- 
basic arc  (p,  q)  with  respect  to  the  arcs  in  the  basis  graph.  The  calculation  can 
be  accelerated  by  tracing  predecessor  paths  to  their  intersection,  or  to  loop  nodes 
(if  the  intersection  occurs  on  a loop),  and  performing  a single  aggregated  calcu- 
lation (rather  than  a sum)  from  the  point  of  intersection.  (If  desired,  one  may 
re-scale  arcs  in  P^  and  P^  that  are  mirror  arcs  to  obtain  the  representation  co- 
efficients in  terms  of  the  original  arcs.  However,  this  is  unnecessary  in  the 
execution  of  the  extreme  point  algorithm.) 


4 . 3 Updating  Maintaining  t hj?  EAPl  Bas  is  St ructure 

To  develop  the  special  rules  and  relationships  underlying  the  complete 

basis  exchange  operation,  subject  to  maintaining  the  EAPI  orientation,  we  will 

make  use  of  the  following  definitions.  From  Remark  4,  the  basic  arcs  with  a 

nonzero  coefficient  in  the  representation  of  a nonbasic  arc  (p,  q)  are  contained 

in  P and  P . These  arcs  can  be  separated  into  three  mutually  exhaustive  sets 

P q 

at  least  one  of  which  is  nonempty: 

a)  PpDPq 

(2)  P - P 

p q 

(3)  P - P 

q p 


If  i is  some  node  in  the  predecessor  path  of  node  j,  we  define  P^^  to  be  the  set 

of  arcs  in  the  predecessor  path  from  node  j to  node  i.  If  i = j or 

i is  not  an  ancestor  of  j,  P^^  = 0.  Then,  by  the  previous  definition  of  0(S) 

( = 1/product  of  the  multipliers  of  the  arcs  in  set  S)),  we  have  “ l/(the 

product  of  the  multipliers  of  the  arcs  in  P^  and  not  in  P^)*  (Recall  by  the 

convention  0(0)  = 1,  that  6(Pj^j)  = 1 if  i = j or  i is  not  an  ancestor  of  j.) 

We  may  further  elaborate  the  properties  of  the  function  0 (to  facilitate 
subsequent  derivations  involving  sets  of  arcs  on  path  segments,  such  as  P^,  P^j > 
etc.)  as  follows. 


Remark  5 : If  , 


, is  any  partition  of  a set  of  arcs  R,  then 


9(R)  = 0(R^)0(R2) 


0(R  ). 

V 


In  particular,  if  k is  a node  on  the  predecessor  path  of  node  j,  and  i is  a node 


on  the  predecessor  path  of  node  k,  so  that  P = P.,IJP,  •>  then  0(P.,)  = 

ij  ik'-'  kj  ij 

0 (Pfk)  6 (Pkj ) • allow  k = i or  k = j , in  which  case  tir  0(P^^)  = 1 

and  the  relation  still  holds.) 


22 


Remark  6 : If  R*  is  the  set  of  mirror  arcs  of  the  arcs  in  R,  then  9(R*)  = 1/0(R). 

(Hence  in  particular  6(P.,*)  = 1/0(P.,). 

ij  ij 

We  will  continue  to  use  the  asterisk  to  indicate  a "mirror  set"  in  the 
following.  Also,  for  completeness,  it  is  useful  to  define  0*  = 0 (hence  6(S)  = 

0(S*)  = 1 if  S = 0).  However,  to  avoid  cluttered  notation  such  as  0({(i,  j)})  ; 

and  0({(i,  j)}*),  we  will  let  (1,  j)  represent  both  a specified  arc  from  i to  j i 

and  the  set  consisting  of  this  arc. 

The  representation  of  an  incoming  arc  (p,  q)  is  by  definition  a weighted  ’ 

linear  combination  of  the  arcs  in  an  EAPI  basis  graph.  By  means  of  this  repre- 
sentation, any  flow  change  in  (p,  q)  induces  a "corresponding"  flow  change  in 

each  arc  of  its  representation.  Thus,  in  particular,  if  A is  the  flow  change  in 
(p,  q),  and  Y is  the  coefficient  of  a given  arc,  then  -Ay  is  the  flow  change  on  the 
basic  arc.  This  means  that  the  flow  change  on  the  basic  arc  will  have  the  same 
sign  as  the  flow  change  on  (p,  q)  if  the  coefficient  of  the  basic  arc  is  negative, 
and  will  have  a flow  change  of  the  opposite  sign  if  the  coefficient  of  the  basic 
arc  is  positive.  An  arc  whose  flow  is  among  those  first  driven  to  an  upper  or 
lower  bound  by  an  attempted  change  of  the  value  of  (p,  q)  (away  from  the  bound 

it  currently  equals)  is  called  a blocking  arc.  (By  our  convention  of  employing  ; 

mirror  arcs,  (p,  q),  being  nonbasic,  always  equals  its  lower  bound  prior  to  the 

attempted  flow  change.  It  should  be  noted  that  (p,  q)  can  be  a blocking  arc.) 

In  the  case  of  a degenerate  pivot,  the  flow  on  a blocking  arc  currently 
equals  one  of  its  bounds  and  would  violate  that  bound  with  any  change  in  the  flow 
of  the  incoming  arc  (p,  q).  To  maintain  primal  feasibility,  the  outgoing  arc 
must  always  be  one  of  the  blocking  arcs. 

We  consolidate  these  and  several  earlier  observations  by  defining: 
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Q = the  sec  of  basic  arcs  with  positive  representation  coefficients  for  (p,  q) 

P = the  set  of  basic  arcs  with  negative  representation  coefficients  for  (p,  q) 
i 


(where,  as  will  be  seen,  Q is  associated  with  node  q and  P is  associated  with  node 
P). 

Remark  7 : Q is  the  set  of  arcs  whose  flows  are  decreased  and  P\J(p,  q)  is  the 

set  of  arcs  whose  flows  are  increased  when  the  flow  of  the  incoming  arc  (p,  q)  is 

Increased.  Further,  PUQcp  fjP 

P q 

Remark  8:  The  outgoing  arc  (r,  s)  satisfies: 

(r,  s)eQ  if  (r,  s)  is  blocking  at  its  lower  bound. 

(r,  s)ePU(p,  q)  if  (r,  s)  is  blocking  at  its  upper  bound. 

We  are  now  in  a position  to  state  the  major  result  of  [18]  which  reduces  a 
diverse  collection  of  basis  updating  configurations  to  two  simple  cases  and 
identifies  the  appropriate  reorientation  required  to  maintain  the  EAPI  basis 
structure . 

Theorem  1 [18]  : For  a basis  exchange  involving  the  incoming  arc  (p,  q)  and  the 

outgoing  arc  (r,  s) , the  EAPI  rooted  loop  orientation  is  maintained  as  follows: 

Reverse  the  orientation  of  the  arcs  of  P U(p,  q)  if  (r,  s)cPU(p,  q) 

sp 

Reverse  the  orientation  of  the  arcs  of  P if  (r,  s)EQ. 

sq 

It  may  be  noted  that  reversing  the  orientation  of  a set  of  arcs  S corresponds 

to  replacing  S by  its  mirror  set  S*.  To  translate  and  extend  the  implications 

of  this  important  result  still  more  fully  to  the  mirror  arc  setting,  and  set  the 

stage  for  further  developments,  we  define  for  the  situation  in  which  P H ^ ^ 

[ P q 

• X = the  unique  node  on  the  predecessor  path  of  node  p such  that 

i P * P - P 

! xp  p q 

I I y = the  unique  node  in  the  predecessor  path  of  node  q such  that 


P = P - P 

yq  q p 
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z = the  first  intersection  node  of  P and  P , i.e.,  the  unique  node  of 

sp  sq 


these  paths  such  that  P OP 

sp'  sq 

and  P (JP  = P IJP  - (P  HP  ) 

zq  zp  sq'^  sp  sq'  ' sp 

Figures  4 and  5 illustrate  these  definitions. 


P , or  equivalently  P Op  = 0 
sz  zp  zq 


Note  that  when  x = y,  node  z = node  x = node  y. 

)eP 


When  X y 


r X if  (r,  s) 
' ^ Lyif(r,  s) 


eP 


yx 

xy 


Remark  9:  If  P H?  ^ 0,  P = P IIP  and  0 = P I lO  for  some  P , Q , possibly 

— — — p q xp^  o yq^  o o o 

empty  and  contained  in  P HP 

P q 

The  significance  of  the  identification  of  node  z is  demonstrated  by  the 
following  consequence  of  Theorem  1 which  refines  the  specification  of  the  sets  P 
and  Q as  follows: 

Corollary  1:  A new  loop  is  created  by  adding  the  incoming  arc  (p,  q)  and  dropping 

the  outgoing  arc  (r,  s)  if  and  only  if  (r,  s)eP  HP  ^nd  the  composition  of  the 

P q 

new  loop  in  the  EAPI  basis  is  given  by 

(P,  if  s)e(P  nP  )0P 

^4  P 4 

(p,  q)LJP2pUP2q*  if  s)e(Ppnp^)nq 

Note  that  the  composition  of  the  loop  when  (r,  s)eP  is  exactly  the  mirror  of  its 
composition  when  (r,  s)£Q. 


5.  THE  NEW  ALGORITHM 

The  preceding  development  consolidates  many  of  the  fundamental  results  for 
GN  problems  that  can  be  found  piecemeal  in  a variety  of  rexerences  [1,  7,  12,  16, 
18,  22,  24].  By  means  of  the  generalized  mirror  arc  concept,  several  useful  inter- 
relationships have  been  made  visible  that  are  not  explicit  in  the  literature.  Our 
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primary  motivation  for  this  development,  however,  is  not  simply  to  bring  together 
I scattered  concepts  in  a unified  format,  but  to  lay  a foundation  for  character- 

I izing  the  special  structural  attributes  and  convergence  properties  of  the  new 

algorithm.  The  following  sections  are  devoted  to  this  end. 

5 . 1 Fundamental  Definitions 

An  arc  has  lower  leeway  if  its  flow  is  strictly  greater  than  its  lower  bound, 
and  has  upper  leeway  if  its  flow  is  strictly  less  than  its  upper  bound.  An  arc 
which  has  both  lower  and  upper  leeway  will  be  called  a double  leeway  arc.  It 
follows  that  an  arc  has  lower  leeway  if  and  only  if  its  mirror  has  upper  leeway, 
and  vice  versa.  Consequently,  if  an  arc  is  a mirror  arc,  then  it  has  lower  or 
upper  leeway  according  to  whether  the  original  arc  has  the  opposite  type  of 
leeway . 

A basis  B is  defined  to  be  a strongly  convergent  (SC)  basis  if  and  only  if 
the  EAPI  basis  graph  satisfies  each  of  the  following  conditions: 

(1)  Every  arc  has  lower  leeway. 

(2)  The  loop  factor  0(L)  of  each  loop  L of  the  basis  graph  other  than  a 
self-loop  satisfies  0(L)<1. 

5 . 2 ^ Basis  Equivalence 

Consider  the  basis  shown  in  Figure  2.  At  least  one  of  the  arcs  with  0 flow 
is  an  original  arc,  thus  this  is  not  an  SC  basis  since  not  all  arcs  have  lower 
leeway.  In  general,  it  is  clear  that  many  bases  for  generalized  networks  will 
not  be  SC  bases. 

A necessary  condition  for  a basis  to  be  equivalent  to  an  SC  basis  is  that  all 
nonloop  arcs  satisfy  Condition  1 of  the  SC  basis  definition.  A sufficient  condi- 
tion for  equivalence  is  given  by  the  following  remark. 


t m ■ L 


.J 
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Remark  10:  A basis  that  satisfies  Condition  1 of  the  SC  basis  definition  can 

be  reoriented  to  create  an  SC  basis  if  and  only  if  the  arcs  of  each  loop  L with 
9(L)>1  are  double  leeway  arcs. 

Proof ; If  0(L)  = 1 and  L is  not  a self-loop,  the  vectors  associated  with  the 
arcs  in  L form  a linearly  dependent  set  as  a direct  consequence  of  the  solution 
identified  in  Remark  3,  which  is  undetermined  when  6(L)  = 1.  Consequently  such 
an  L cannot  be  in  the  basis.  Thus,  we  need  only  consider  the  case  where  6(L)>1. 
When  L is  oriented  in  the  reverse  direction,  L is  replaced  by  its  mirror  set  L* 
and  6(L*)  = (1/9(L))<1  by  Remark  6.  Since  the  mirror  of  an  arc  has  lower  leeway 
if  and  only  if  the  arc  has  upper  leeway,  the  conclusion  follows  at  once. 

By  the  same  reasoning,  if  the  only  arcs  without  lower  leeway  are  loop  arcs, 
the  basis  can  be  reoriented  to  create  an  SC  basis  if  and  only  if  every  basis  loop 
L that  contains  such  an  arc  satisfies  0(L)>1  and  consists  entirely  of  arcs  with 
upper  leeway. 

It  should  be  pointed  out,  however,  that  it  is  always  possible  to  find  an 
initial  feasible  SC  basis  (which  may  involve  artificial  vectors  for  a GN  pro- 
blem). The  following  remark  establishes  this  statement. 

Remark  11 : A basis  that  consists  of  positive  and  negative  unit  vectors  and 

-E^,  where  E^  represents  the  i*”^  column  of  the  identity  matrix,  is  an  SC  basis 

provided  E^  is  basic  if  b . >0  and  -E^  is  basic  if  b.  <0  (where  b.  is  the  i*’^ 

1 1-1 

component  of  the  right-hand  vector  b of  (3)). 

Proof : The  proof  follows  immediately  from  the  definition  of  an  SC  basis  and  the 

comments  regarding  the  EAPI  orientation  of  surplus  arcs. 

5.3  Special  Properties  of  an  SC  Basis  and  the  SC  Algorithm 


In  this  section,  we  will  show  that  every  SC  basis  exhibits  two  fundamental 
properties : 
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(1)  The  signs  of  the  nonzero  coefficients  in  the  representation  of  a 
nonbasic  arc  with  respect  to  an  SC  basis  can  be  determined  prior  to  actually 
finding  the  representation  (by  means  of  a complete  characterization  of  P and 


(2)  There  exists  a unique  arc  in  the  EAPI  basis  graph  which  must  be  chosen 
to  leave  the  basis  at  each  basis  exchange  step  if  the  SC  basis  structure  is  to 
be  maintained. 

These  properties  have  important  implications  for  the  computational  efficiency 
of  the  new  algorithm  founded  on  the  SC  bases  (hereafter  called  the  SC  algorithm). 

A further  useful  property  that  is  immediately  available  from  the  definition 
of  an  SC  basis  and  the  preceding  results  follows. 

Remark  12 : A basis  exchange  step  executed  relative  to  an  SC  basis  will  be  non- 

degenerate if  the  outgoing  arc  (r,  s)£(p,  q)UQ. 

We  now  characterize  the  composition  of  P and  Q for  an  SC  basis. 

Theorem  2:  In  an  SC  basis,  ilPplP  =0,  P=P  and  Q = P . Otherwise,  P 

P q P q 

consists  of  consecutive  arcs  in  the  predecessor  path  of  node  p and  Q consists  of 

consecutive  arcs  in  the  predecessor  path  of  node  q,  such  that 

P = P JP 
xp^  o 

Q = P UQ 

yq  o 


P = 0 or  P or  P np 

o yx  p' ' q 

Q = 0orP  orPflP 
o xy  P q 
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Further,  for  x y 


P C 

xy 


P if  and  only  if  a 6(P  )<6(P  ) 

o pq  yq  yp 

Q if  and  only  if  a 0(P  )>6(P  ) 

' o pq  yq  yp 


P C 

yx  ^ 


P if  and  only  if  a 6(P  )<6(P  ) 
o pq  xq  xp 

Q if  and  only  if  a 6(P  )>6(P  ) 
o pq  xq  xp 


and  for  x = y 


P if  and  only  if  a 9(P  )<6(P  ) 

P Qp  o pq  yq  yp 

Q if  and  only  if  a 9(P  )>9(P  ) 
o pq  yq  yp 


Proof : The  essential  content  of  the  theorem  not  previously  established,  is  the 

identification  of  the  form  of  P and  Q for  the  SC  basis  when  P PlP  ^ 0-  Let  L 

o o P q 

denote  the  unique  loop  contained  in  P HP  • Allowing  either  or  both  of  P and 

P q xp 

P to  be  empty  (i.e.,  for  x = p or  y = q , in  which  case  9(P  ) or  9(P  ) = 1), 

yq  I'  :>  'ir  xp  yq 

the  representation  coefficient  y for  any  arc  (i,  j)eP^PlP^)  can  be  characterized 

from  Remarks  3 and  4,  as  follows: 

Y = a 9(P.  ) - 9(P.  ) if  (i,  i)^L  (which  can  occur  only  if  x = y)  or  if 
(i,  j)  is  the  slack  arc  of  a self-loop. 

Y = ^^^pq®^^iq^  ~ 9(P^p))/(l  - 9(L))  if  (i,  j)e  L and  L is  not  a self-loop. 
Because  0(L)<1  for  the  SC  basis,  y can  be  written  y=a  (a  9(P.  )-0(P.  )) 

pq  iq  ip 

where  a>0,  in  both  instances.  If  x = y,  the  results  follow  Immediately  since 

(1)  P = P =0  and  (2)  the  sign  of  the  coefficients  on  all  the  arcs  in 

PpOPq  correspond  to  the  sign  of  the  transmitted  requirement  to  node  x = y 

which  is  a 9(P  )-9(P  ) = a 9(P  )-9(P  ).  If  x y,  first  suppose  (i,  j)eP 

pq  xq  xp  pq  yq  yp  xy 


^4 


I- 


r 
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The 


n we  can  write  P.  = P,  I JP  and  P.  = T.  IIP  . Hence 
1C|  ly^  yq  ip  ly'^  yp 


Y = a0(P.  ) |a  0(P  )-0(P  ) 

ly  [_pq  yq  yp 


] 


from  which  it  follows  that  the  sign  of  Y depends  only  on  the  sign  of  the 

quantity  in  braces,  and  is  independent  of  the  choice  of  arc  (i,  j)  itself,  as 

long  as  (1,  i)eP  . Hence  P r"  P or  P Cl  Q according  to  the  inequality 
xy  xy^  o xy^  o 

stipulated  in  the  theorem. 

Next,  if  (i,  i)eP  , we  can  write  P.  = P.  MP  and  P.  = P.  IIP 
’ ' > yx’  iq  xq  ip  ix^  xp 

Consequently , 


Y = CX0(P.  ) 

IX 


a 0(P  )-0(P  ) 

pq  xq  xp 


and  by  corresponding  reasoning  we  conclude  or  Pyj^  ^y  Inequalities 


specified  in  the  theorem. 


Finally,  it  is  necessary  to  show  that  P^  and  (hence  P and  Q)  are  indeed 
sets  of  consecutive  arcs  in  predecessor  paths,  i.e.. 


P = 

0 or 

P 

or  p n P 

and  Q = 1 

0 

yx 

p q 

o 

This 

assertion 

is  equivalent 

to  showing 

P P 

0 xy 

the 

foregoing 

results 

and 

the 

fact  that 

'■on'Jo  ■ 

the 

expression 

for  Y 

when 

(1, 

j)eP  by 

yx 

yx 

0(P 

)0(P 

) 

= 0(p  )0(p  )0(p 

yx 

xq 

yx 

xy  yq 

0(P 

)0(P  ) = 

0(P  ^),  ' 

yx  xp 

yp 

0(P 

)Y 

= a0(P.  ) 

a 0(P 

yx 

IX 

pq  yq 

xy  P q 

Q ^ P (on  t 
o yx 


yq 


Thus  the  sign  of  y for  (i,  j)eP  depends  on  the  quantity  in  the  braces  in  this 

latter  expression.  Comparing  this  quantity  to  the  quantity  in  the  braces  in  the 

expression  for  Y when  (i,  J)gP  , and  noting  0(L)<],  it  follows  that  y>0  for 

xy  - 

(i,  j)eP  implies  Y^O  for  (1,  i)eP  , and  this  together  with  its  contrapositive 
yx  xy 
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leads  to  the  conclusion  of  the  theorem. 

The  preceding  characterizations  make  it  possible  to  state  the  SC  algorithm 
and  prove  its  special  properties.  For  this  purpose,  we  suppose  the  arcs  of  P 
and  Q are  indexed  in  ascending  order,  in  the  same  sequence  as  they  are  encountered 
by  a trace  of  the  predecessor  indexing,  starting  at  nodes  p and  q.  (Thus,  if  P is 
not  empty,  (p(p),  p)  is  the  first  arc  of  P,  and  if  Q is  not  empty,  (p(q),  q)  is  the 
first  arc  of  Q.)  In  addition,  we  will  augment  the  set  P by  the  addition  of  arc 
(p,  q)  , to  give  the  set  (p,  q)UP,  and  treat  arc  (p,  q)  as  the  first  arc  of  the 
augmented  set  (whereupon  the  first  arc  of  P is  the  second  arc  of  this  set,  etc.). 

The  complete  form  of  the  SC  algorithm  may  be  described  as  follows.  (The 
specifications  for  determining  and  updating  node  potentials  in  this  description 
can  be  implemented  efficiently  by  reference  to  the  procedures  described  in  [16, 
18].) 

The  SC  Algorithm 

1.  Begin  with  an  SC  Basis  (Remark  11)  and  determine  node  potentials  for  each 

node  k,  so  that  reduced  cost  a . .TT  . -it  . -c . . = 0 for  each  basic  arc  (1,  i). 

11  1 1 ij 

2.  Select  the  incoming  arc  (p,  q)  to  be  any  nonbasic  arc  whose  reduced  cost 

a TT  -TT  -c  is  profitable  (i.e.,  positive).  If  no  such  arc  exists  (and  all 
pq  q p pq 

artificial  arcs  have  zero  flow),  the  problem  is  solved  and  the  current  arc  flows 
are  optimal. 

3.  Select  the  outgoing  arc  (r,  s)  to  be; 

(1)  the  last  blocking  arc  in  Q if  any  blocking  arc  is  a member  of  that  set. 

(2)  the  first  blocking  arc  of  (p,  q)(jP  if  there  is  no  blocking  arc  in  Q. 
h.  Execute  a basis  exchange  step,  reorienting  the  basis  by  Theorem  1 and 
determining  the  updated  flows  and  node  potentials.  Then  return  to  instruction  2. 


*■  ■ 
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Recall  that  the  incoming  and  outgoing  arcs  can  be  self-loops  (i.e.,  (r,  s)  can 
be  a slack  arc  and  (p,  q)  can  be  either  a slack  or  a surplus  arc).  By  convention, 
for  these  cases,  we  define  7T.^^  = 0 and  = 0 (for  slack  and  surplus  arcs, 

respectively).  This  yields  correct  reduced  cost  expressions  for  self-loops. 
Theorem  3:  The  SC  algorithm  maintains  the  SC  basis  structure  at  each  basis  ex- 

change step.  Moreover,  any  other  choice  of  (r,  s)  destroys  this  structure. 

Proof:  First,  we  show  that  all  arcs  of  (p,  q)e  PUQ  (the  only  arcs  whose  flows 
change)  will  have  lower  leeway  after  the  exchange  step  if  and  only  if  (r,  s)  is 
selected  from  this  set  as  specified.  If  (r,  s)eQ,  the  flow  change  is  nonde- 
generate by  Remark  12.  None  of  the  arcs  in  Q indexed  higher  than  (r,  s)  are 
blocking  arcs.  Hence  these  arcs  retain  their  lower  leeway  and  further  are  not 
reoriented  by  the  basis  exchange  step.  If  a blocking  arc  existed  that  did  have  a 
higher  index  than  (r,  s),  it  would  not  have  lower  leeway  after  the  pivot  and 
would  also  not  be  reoriented,  thereby  violating  the  SC  structure.  The  arcs  of 

which  are  reoriented  (Theorem  1)  must  have  upper  leeway  before  reorientation 

due  to  the  fact  that  P Cl  Q and  the  flow  change  is  nondegenerate.  Consequently, 

sq 

after  reorientation  they  have  lower  leeway.  Finally,  arcs  in  (p,  q)e  P,  which 
are  not  reoriented,  have  lower  leeway  because  their  flows  are  increased.  If 
(r,  s)e(p,  q)lJP,  then  none  of  the  arcs  in  (p,  q)lJP  indexed  lower  than  (r,  s) 
are  blocking,  and  hence  must  have  upper  leeway.  These  are  precisely  the  arcs 
that  are  reoriented  (Theorem  1)  and  hence  have  lower  leeway  after  reorientation. 

If  any  arc  indexed  lower  than  (r,  s)  had  been  blocking,  then  after  reorientation 
such  an  arc  would  have  had  no  lower  leeway,  thereby  violating  the  SC  structure. 

No  arcs  in  Q are  blocking,  and  hence  retain  lower  leeway,  and  arcs  in  P not  re- 
oriented have  lower  leeway  since  their  flows  are  only  increased  or  remain  the 


same . 


Now,  we  must  establisti  that  0(L)<1  for  any  loop  L created  by  the  basis  ex- 


lBttlLS.S](£:£2jBi&iGS£^9[iCDd£llAX 
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change  step,  other  than  a sell-Joop.  By  Corollary  1,  we  need  only  consider  the 

case  where  (r,  s)£(P  flP  )/!  P and  (r,  s)e(1>  f|!’  )PlQ-  Prom  Section  4.3,  if 

p q p ' q ' ' ' 

X = y,  z = X = y.  If  X ^ y,  z = X if  (r,  s)cP  and  z = y if  (r,  s)eP 

yx  xy 

Assuming  that  (r,  s)eP,  tlie  new  loop  L formed  is  (p,  q)*ljP  .I  I P 

zp*^  zq 

If  X = y,  by  Theorem  2 and  the  substitution  of  z,  a 0(P  )<0(P  ) which  implies 

pq  zq  zp 

9(L)<1.  If  X y,  by  Theorem  2,  we  need  only  consider  the  case  (r,  s)eP  . By 

yx 

Theorem  2 and  the  substitution  of  z,  a 0(P  )<9(P  ) which  again,  implies 

pq  zq  zp 

0(L)<1.  Similar  reasoning  applies  when  (r,  s)eQ. 

5 . 3 Convergence  Propert ies  of  the  SC  Algorithm 

Having  established  the  way  to  create  and  maintain  an  SC  basis,  and  having 
characterized  the  sufficiency  condition  by  which  a pivot  in  this  basis  is  non- 
degenerate, we  will  now  show  that  the  SC  algorithm  is  finitely  convergent  without 
any  reliance  upon  "external"  techniques  such  as  perturbation.  Moreover,  we  will 
show  that  the  convergence  through  a succession  of  degenerate  pivots  is  of  a much 
stronger  form  than  that  established  by  reference  to  "perturbation"  by  showing 
that  the  dual  variables  (node  potentials)  which  change  after  a basis  exchange 
step  change  in  a uniform  direction  throughout  any  sequence  of  degenerate  pivots. 

Removing  the  outgoing  arc  (r,  s)  before  adding  the  incoming  arc  (p,  q) 
always  creates  a unique  tree  in  the  EAPI  basis  graph.  (This  tree  will  contain  all 
nodes  of  a former  quasi-tree  if  the  outgoing  arc  is  a loop  arc.  The  tree  may 
also  consist  of  only  a single  node.)  We  let  T denote  the  set  of  nodes  in  this 
tree.  The  only  node  potentials  which  change  in  the  new  basis  are  those  associated 
with  the  nodes  in  T.  A change  in  the  node  potential  of  any  particular  node  of  T 
induces  a change  in  the  node  potentials  of  all  its  successor  nodes.  This  change 
is  characterized  as  follows. 

Remark  13:  In  order  to  maintain  complementary  slackness  (reduced  arc  costs  of  0) , 

a change  of  a in  the  node  potential  of  node  i in  T induces  a change  of  aO(Pj.)  in 
the  node  potential  of  any  successor  node  j in  T. 
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Proof : Assume  tliat  the  nodes  in  the  predecessor  path  from  node  j to  node  i 

are  numbered  beginning  at  node  i as  i,  i+1 , . . . , j.  Let  Tij^  denote  the  current 
node  potential  and  denote  the  change  in  this  node  potential  for  node  k.  In 
order  for  complementary  slackness  to  hold, 

where  c^^  denotes  the  cost  attached  to  the  arc  (k,k+l),  k = i,  . . . , j-l. 

The  current  node  potentials  also  satisfy  complementary  slackness: 

-\+'>k+i\+r'k,k+i  ■ “ 


for  each  arc  (k,k+l),  k = i. 


, j-l- 


Combining  (11)  and  (12)  we  obtain 


The  remark  follows  at  once. 

Remark  14 : A change  in  the  node  potential  of  node  i induces  a change  of  the 

same  sign  in  the  node  potentials  of  all  its  successors  in  T if  the  conditions  of 
complementary  slackness  are  maintained. 

Proof : A direct  consequence  of  Remark  13. 

Using  the  above  remarks,  we  will  now  prove  that  the  SC  algorithm  is  finitely 
Convergent  and  this  convergence  has  a particularly  strong  character.  Let  it 
denote  the  vector  of  node  potentials.  An  extreme  point  method  is  said  to  be 
uniformly  converging  if  every  component  of  tt  that  changes  must  strictly  increase 
as  a result  of  a degenerate  pivot.  Clearly  this  form  of  "uniform"  convergence  is 
substantially  stronger  than  the  usual  lexicographic  convergence"  in  which  only 
the  least  indexed  component  of  tt  which  changes  must  increase. 
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Theorem  4:  Every  pivot  of  the  SC  algorithm  causes  the  potentials  of  the  nodes  in  T 

to  change  in  the  following  ways: 

(J)  If  (r,  s)gQ,  the  potentials  of  the  nodes  in  T decrease. 

(2)  If  (r,  s)eP,  the  potentials  of  the  nodes  in  T increase. 

Further,  the  SC  algorithm  is  unifonnly  converging  through  every  sequence  of  de- 
generate pivots. 

Proof : The  incoming  arc  (p,  q)  is  eligible  to  enter  the  basis  if  its  reduced  cost 

is  positive,  i.e. , 


-TT  + a TT 

p pq  q 


c = a >0 

pq 


(13) 


First  suppose  (r,  s)t^P  HP  • By  Theorem  1,  the  tree  associated  with  T will  be 

P q 

rooted  after  re-orientation  at  node  q (node  p)  if  the  outgoing  arc  belongs  to 

P - P (P  - P ).  If  node  q becomes  the  root  (hence  (r,  s)eQ  by  Theorem  1)  the 
q p p q 

potential  of  q must  change  by  6 to  maintain  complementary  slackness,  where 


-IT  + a (it  + 5)-c  = 0 

p pq  q pq 


(14) 


Rewriting  (13)  using  (14)  we  obtain 


r -ct 

6 = < 0. 

a 

pq 

Hence  by  Remark  14,  the  potentials  of  all  nodes  in  T decrease.  If  node  p becomes 
the  root  (hence  (r,  s)eP),  similar  reasoning  discloses  the  potentials  of  nodes  in 
T increase. 

Next  suppose  (r,  s)cP^P|  P^  . Both  node  p and  node  q are  then  in  T and  thus  both 

node  potentials  will  change.  Let  6 represent  the  change  in  ti  and  A represent  the 

P 

change  in  . To  maintain  complementary  slackness  in  the  updated  basis: 


-(tt  +A)+a  (tt  +A)-c  = 0 

p pq  q pq 


(15) 
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It  (r,  s)eQ,  Che  arcs  in  P are  reoriented  in  the  updated  basis  such  that  node  p 

sq 

is  a successor  of  node  q (before  adding  arc  (p,  q),  which  also  makes  q a successor 
of  p).  Thus  from  Remark  13,  a change  of  A at  node  q transmits  a change  of 
A0(P^p)  to  node  p for  the  predecessor  path  that  connects  node  p and  node  q after 
reorienting.  However,  from  (13),  we  obtain  6 = a+a  A which  implies  a+a  A = 

pq  pq 

Ae(p  ). 
qp 


Solving  for  A, 


A = a/(0(P  ) - a ) 

qp  pq 


(16) 


By  Corollary  1,  the  composition  of  P in  terms  of  the  basis  arcs  before  re- 

qp 

orienting  is  ^^pU ^zq*  that  (16)  becomes  upon  rearranging  terms 

A = a 0(P  )/0(P  )-a  0(P  ) . 

zq  zp  pq  zq 

Using  the  appropriate  substitution  for  z when  (r,  s)eP  OP  , x = y;  (r,  s)eP  ; 

p*  ' q xy 

and  (r,  s)sP  , we  obtain  from  Theorem  2 A<0. 
yx 

Since  node  q is  a loop  node  in  the  updated  basis,  all  nodes  whose  potentials 
change  are  successors  of  node  q.  Thus  by  Remark  14,  all  node  potentials  which 
change  will  decrease. 

If  (r,  s)£P, 

6 = - ■ 


1-a  0(P  ) 

pq  pq 

Using  the  same  reasoning  as  above,  it  follows  that  6>0  when  (r,  s)eP  OP  , x = v; 

P ‘ q 

when  (r,  s)eP  ; and  when  (r,  s)eP  and  thus  all  node  potentials  which  change  will 
xy  yx 

increase.  Finally,  since  in  a degenerate  pivot  (r,  s)eP,  the  uniformly  converging 
property  is  established  and  the  theorem  is  proved. 
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6.0  Computat lonal  Advantages 

In  the  SC  algorithm,  the  number  of  possible  degenerate  pivots  is  reduced 
since  the  special  structure  of  the  basis  eliminates  many  of  the  alternative 
basis  representations  for  extreme  paints  normally  examined  by  other  algorithms. 

This  combined  with  the  strong  convergence  properties  of  the  algorithm  should  re- 
sult in  fewer  total  pivots.  Further,  since  no  computational  testing  is  re- 
quired to  ensure  that  only  bases  with  the  SC  structure  are  examined,  there  will 
be  no  accompanying  increase  in  execution  time  per  pivot. 

Hultz,  et  al  [21]  showed  that  by  incorporating  degeneracy  checks  into  a 
GN  code,  solution  times  were  substantially  decreased.  In  this  code,  degeneracy 
checks  are  applied  to  each  arc  traversed.  When  executing  the  SC  algorithm, 
degeneracy  can  occur  only  in  a portion  of  the  basis  graph  (i.e.,  the  set  P) . 

Thus,  degeneracy  checks  need  only  be  applied  to  the  arcs  in  P.  Further  by  the 
pivot  rules  of  the  algorithm,  the  first  arc  encountered  in  P which  will  cause  a 
degenerate  pivot  is  also  the  outgoing  arc;  thus  no  further  tree  traversal  or 
additional  degeneracy  checks  are  required  once  degeneracy  is  detected. 

For  these  reasons,  the  SC  algorithm  when  Implemented  should  result  in 
several  computational  advantages  over  currently  available  codes  [19,  22,  24] 
for  solving  GN  problems. 

7 . 0 Limits  of  Generality 

As  noted  earlier,  methods  with  the  uniform  convergence  property  have  pre- 
viously been  developed  in  the  pure  network  setting  for  assignment,  transportation 
and  transshipment  problems  [2,  3,  4,  11].  There  is  a major  leap  from  the  generality 
of  the  pure  network  setting  to  that  of  the  generalized  network  setting  which,  suc- 
cessfully accomplished  by  the  preceding  development,  raises  the  question  of 
whether  it  is  possible  to  discover  a primal  extreme  point  algorithm  with  the 
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uniform  convergence  property  for  the  still  more  general  case  where  both  ends  of 
a generalized  arc  may  have  multipliers  (i.e.,  nonzero  entries  in  the  associated 
column  of  A)  of  the  same  sign.  In  fact,  our  results  make  it  possible  to  con- 
clusively answer  this  question  in  the  negative,  establishing  the  precise  limits 
to  the  level  of  problem  generality  for  which  such  a method  exists,  and  identify- 
ing the  SC  algorithm  as  the  method  that  attains  that  level  (relative  to  linear 
programming  problems  with  at  most  two  entries  in  each  column  of  A). 

To  demonstrate  this,  consider  first  of  all  the  required  change  in  the  defini- 
tion of  0(S)  for  a set  of  arcs  S when  some  arcs  may  have  multipliers  of  the  same 
sign.  Retracing  the  proof  of  Remark  3,  which  motivates  this  definition  (to  give 
the  form  of  the  representation  coefficients),  it  is  clear  that  S must  specifically 
refer  to  a set  of  arcs  with  an  implied  orientation,  as  induced  by  the  predecessor 


indexing.  Hence,  if  S consists  of  k arcs,  with  multipliers  a » u for  j = 1,  . . . 

i j 

k,  where  a.  is  the  multiplier  at  the  successor  node  and  u,  is  the  multiplier  at  the 
3 J 

predecessor  node  under  the  implied  orientation  (hence  Uj  occupies  the  role  of  the 

-1  multiplier  in  the  preceding  development)  then  the  appropriate  definition  is 
Tl  if  S = 0 
0(S)  = V 1/a^  if  k = 1 


^(1/a  ^)  (-u^/a^  ) (-u^/a^)  . . . (-Uj^_^/aj^  ) if  k > 1 


In  addition,  however,  the  loop  factor  must  be  redefined  to  be  0(L)  instead  of 


0 (L) , where 


0(L)  = (-U  /a  ) . . . (-U  /a  ) 
11  k k 


for  L composed  of  the  arcs  specified  above  to  be  in  S (hence,  in  general,  0(S)  = 

-u  . . -u  0(.S)).  By  these  conventions,  replacing  6(L)  by  0(L)  , the  statement 

1 R—  1 

of  Remark  3 is  valid  for  arcs  with  arbitrary  multipliers. 
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Note  that  when  all  of  the  are  ~1,  as  earlier,  then  Q(S)  = 6(S).  Also, 
for  a self-loop  L,  0(L)  = 0(L)  = 1 by  these  definitions.  Finally,  by  an  obvious 
adaptation  of  the  mirror  arc  concept  to  arcs  whose  multipliers  have  the  same 
sign,  it  can  be  assured  that  = -1  for  all  basic  arcs,  and  the  previous  ex- 
pressions recover  their  validity. 

The  chief  consequence  of  these  observations  is  that  representation  coefficients 
for  the  unit  vector  may  have  both  positive  and  negative  signs  on  the  arcs  of 
the  predecessor  path  of  node  i.  Further,  the  representation  coefficient  of  a 
particular  arc  on  this  path  can  change  its  sign  simply  by  shifting  the  location  of 
node  i on  the  path.  Thus,  the  composition  of  P and  Q,  whose  special  form  is 
essential  to  Theorems  3 and  4,  can  no  longer  be  characterized  in  relation  to  the 
sets  Pp  and  P^  . Finally,  after  a basis  exchange,  the  change  in  the  node  poten- 
tial of  node  j induced  by  a change  of  cx  in  the  potential  of  an  ancestor,  node  i, 
becomes  ot  (Remark  13),  which  may  or  may  not  have  the  same  sign  as  oi  , and 

consequently  it  is  impossible  for  the  node  potentials  to  change  in  any  uniformly 
specifiable  manner  as  required  by  Theorem  4. 

Still  more  simply  (though  at  a less  rigorous  level)  it  is  possible  to  argue 
that  a single  arc  whose  multipliers  are  of  the  same  sign  can  prevent  the  uniform 
convergence  property  from  being  established.  If  such  an  arc  is  nonbasic  and  be- 
comes the  incoming  arc  (p,  q),  then  both  P and  P belong  either  to  P or  Q,  and 

P q 

the  proof  of  Theorem  3,  which  establishes  the  uniqueness  of  the  outgoing  arc 

(r,  s),  discloses  that  there  may  in  this  case  be  irreconcilable  competing  choices 

for  (r,  s)  in  both  P - P and  P - P . Regardless  of  the  choice  of  (r,  s) , 

P q q P 

the  lower  leeway  requirement  is  no  longer  satisfied  by  all  basic  arcs,  and  a sub- 
sequent step  can  therefore  yield  a degenerate  pivot  when  (r,  s)eQ,  invalidating 
the  requirement  needed  to  establish  Theorem  4. 
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:i3  abstract' 
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A major  computational  problem  that  arises  in  the  attempt  to  solve 
generalized  network  and  network-related  problems  is  degeneracy.  In  fact,  using 
primal  extreme  point  solution  techniques,  the  number  of  degenerate  pivots 
performed  frequently  ranges  as  high  as  90%  in  large-scale  applications.  The 
purpose  of  this  paper  is  to  develop  a special  set  of  structural  and  logical 
relationships  for  generalized  network  problems  that  yields  a new  primal  extreme 
point  algorithm  which  avoids  and/or  exploits  degeneracy.  The  major 
mathematical  differences  between  this  new  algorithm  and  the  simolex  method 
are  (1)  each  basis  examined  is  restricted  to  have  a special  topology.  (2)  the 
algorithm  is  finitely  convergent  without  reliance  upon^external^^echniques  such 
as  lexicography  or  perturbation,  and  (3)  a special  screening  criterion  for 
nondegenerate  basis  exchanges  is  available  that  allows  some  of  these  exchanges 
to  be  recognized  prior  to  finding  the  representation  of  an  incoming  arc.  For 
these  reasons,  this  algorithm  has  several  computational  advantages  over  the 
simplex-based  codes  recently  developed  for  solving  generalized  network 
problems. 
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